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Using a sample of 9.7 x 10 BB meson pairs collected with the CLEO 

pq \ detector, we study B decays to the Xci and Xc2 charmonia states, which 

are reconstructed via their radiative decays to J/tp. We first measure the 
branching fraction for inclusive Xci production in B decays to be B{B — > 

^|- . XciX) = (4.14 ± 0.31 ± 0.40) x 10~ 3 , where the first uncertainty is statisti- 

cal and the second one is systematic. We derive the branching fractions for 

Q^ | direct Xci and Xc2 production in B decays by subtracting the known contri- 

bution of the decay chain B — » ip(2S)X with i/j(2S) — » Xcl,27- We obtain 
B[B -> Xd(direct)X] = (3.83 ±0.31 ±0.40) x 10~ 3 . No statistically significant 
signal for Xc2 production is observed in either case. Using the Feldman-Cousins 
approach, we determine the 95% confidence intervals to be [0.2, 2.0] x 10 -3 
for B(B -> Xc2X), [0.0,1.7] x 10~ 3 for B[B -» Xc2 (direct )X], and [0.00,0.44] 
for the ratio T[B — > Xc2(direct)X] /L[i? — ► Xci(direct)X]. We also measure 
the branching ratio T[B — > x c2 (direct)X s ]/TL9 — ► x c i(direct)X,] for different 
X s configurations by reconstructing B decays into exclusive final states with 
J/ip, 7, a kaon, and up to four pions. For all the X s configurations we ob- 
serve a strong Xci signal yet no statistically significant Xc2 signal. We discuss 
how our results compare with theoretical predictions in different models of 
charmonium production. 
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The recent measurements of charmonium production in various high-energy physics reac- 
tions have brought welcome surprises and challenged our understanding both of heavy-quark 
production and of quarkonium bound state formation. The CDF and DO measurements ]I[ 
of a large production rate for charmonium at high transverse momenta (Pt) were in sharp 
disagreement with the then-standard color-singlet model. The development of the NRQCD 
factorization framework has put the calculations of the inclusive charmonium production 
on a rigorous footing. The high-P^ charmonium production rate at the Tevatron is now well 
understood in this formalism. The recent CDF measurement of charmonium polarization [[|, 
however, appears to disagree with the NRQCD prediction. The older color-evaporation 
model accommodates both the high-P^ charmonium production rate and polarization mea- 
surements at the Tevatron Q. 

Inclusive B decays to charmonia offer another means by which theoretical predictions 
may be confronted with experimental data. The color-singlet contribution, for example, is 
thought to be § a factor of 5-10 below the observed inclusive J/ip production rate |J. A 
measurement of the Xc2-to-x c i production ratio in B decays provides an especially clean test 
of charmonium production models. The V — A current 07^(1 — 75)0 cannot create a cc pair in 
a 2S+1 Lj = 3 P2 state, therefore the decay B —>■ Xc2X is forbidden at leading order in a s in the 
color-singlet model . The importance of the color-octet mechanism for Xc production in B 
decays was recognized || even before the development of the NRQCD framework |2|]. While 
the NRQCD calculations cannot yet produce sharp quantitative predictions for the Xc2-to- 
Xd production ratio in B decays ||, we can consider two limiting cases. If the color-octet 
mechanism dominates in B — > XcjX decays, then the Xc2-to-Xci production ratio should be 
5:3 because the color-octet contribution is proportional to 2 J + 1. In contrast, if the color- 
singlet contribution dominates, then Xc2 production should be strongly suppressed relative 
to Xci production. The color-evaporation model predicts the ratio to be 5 : 3 [f|. 

Our data were collected at the Cornell Electron Storage Ring (CESR) with two configu- 
rations of the CLEO detector called CLEO II [IIJ and CLEO II. V 0. The components of 
the CLEO detector most relevant to this analysis are the charged particle tracking system, 
the Csl electromagnetic calorimeter, the time-of-flight system, and the muon chambers. In 
CLEO II the momenta of charged particles are measured in a tracking system consisting of 
a 6-layer straw tube chamber, a 10-layer precision drift chamber, and a 51-layer main drift 
chamber, all operating inside a 1.5 T solenoidal magnet. The main drift chamber also pro- 
vides a measurement of the specific ionization, dE/dx, used for particle identification. For 
CLEO II. V, the straw tube chamber was replaced with a 3-layer silicon vertex detector, and 
the gas in the main drift chamber was changed from an argon-ethane to a helium-propane 
mixture. The muon chambers consist of proportional counters placed at increasing depths 
in the steel absorber. 

We use 9.2 fb _1 of e + e~ data taken at the T(45) resonance and 4.6 ftr 1 taken 60 MeV 
below the T(4S) resonance (off-T(45 l ) sample). Two thirds of the data were collected with 
the CLEO II. V detector. The simulated event samples used in this analysis were generated 
with a GEANT-based [O] simulation of the CLEO detector response and were processed in 



a manner similar to the data. 

We reconstruct the Xci,2 radiative decays to Jfij). The branching fractions for the Xci,2 — > 
J/ip 7 decays are, respectively, (27.3±1.6)% and (13.5±1.1)%, whereas the branching fraction 
for the Xco — ► J/4>1 decay is only (0.66 ± 0.18)% flj~4] . In addition, the Xco production rate 



in B decays is expected to be smaller than the Xci,2 rates |5||§. We therefore do not attempt 
to measure Xco production in this analysis. 

The J/ip reconstruction procedure is described in Ref. [[Uj and summarized here. We 
reconstruct both J /if) — > /i + \i~ and J /if) — > e + e~ decays, recovering the bremsstrahlung 
photons for the J /if) — > e + e~ mode. We use the normalized invariant mass for the J /if) 
candidate selection (Fig.l of Ref. |13|| ). For example, the normalized J /if) — ► fi + fi~ mass 
is defined as [M(/x + /x~) — Mj/^]/a(M), where Mj/^ is the world average value of the J /if) 
mass [14[] and a(M) is the expected mass resolution for that particular /i + yU~ combination 
calculated from track four-momentum covariance matrices. We require the normalized mass 
to be between —6 and +3 for the J/ip —> e + e~ candidates and between —4 and +3 for the 
J /if) — > n + fi~ candidates. The momentum of the J/ip candidates is required to be less than 
2 GeV/c, which is slightly above the maximal J/ip momentum in B decays. 

Photon candidates for x&,2 - > J/ipl reconstruction must be detected in the central 
angular region of the calorimeter (| cos# 7 | < 0.71), where our detector has the best energy 
resolution. Most of the photons in T(AS) — »■ BB events come from tt° decays. We therefore 
discard those photon candidates which, when paired with another 7 in the event, produce a 
normalized ir° —>■ 77 mass between —3 and +2. 

In the first part of this work, called the inclusive analysis, we investigate B — > Xci^X 
decays reconstructing only J/ip and 7. We determine the xa and Xc2 yields in a binned 
maximum-likelihood fit to the mass-difference distribution M^J/tpj) — M(J/ip) (Fig. [l|a), 
where M(J/ip) is the measured mass of a J/ip candidate. The excellent electromagnetic 
calorimeter allows us to resolve the xa and Xc2 peaks. The M(J/ipj) — M(J/ip) mass- 
difference resolution is 8 MeV/c 2 and is dominated by the photon energy resolution. The 
bin width in the fit is 1 MeV/c 2 . The background in the fit is approximated by a 5th-order 
Chebyshev polynomial, chosen as the minimal-order polynomial well fitting the background 
in a high-statistics sample of simulated T(45) — > BB events. All the polynomial coefficients 
are allowed to float in the fit. The xa and Xc2 signal shapes are fit with templates extracted 
from Monte Carlo simulation; only the template normalizations are free in the fit. The Xa 
and Xc2 signal yields in the T(4S r ) data are iV ON (x c i) = 672 ± 47(stat) and iV ON (x c2 ) = 
83 ± 37(stat). The xa and Xc2 yields in off-T(4S) data are both consistent with zero: 
N OFF (xd) — 4 ± 7(stat) and iV OFF (x C 2) = 1 ± 7(stat). Subtracting the contributions from 
non-BB continuum events, we obtain the total inclusive B — > XaX and B — > Xc2X event 
yields N(B -> XaX) = 664 ± 49(stat) and N(B -> X c2X) = 81 ± 39(stat). 

Taking into account the systematic uncertainties associated with the fit, we determine 
the B — > Xc2X signal yield significance to be 2.0 standard deviations (a). Subtracting the 
known contribution of the decay chain B — > ip(2S)X with ip(2S) —>■ Xc2l and accounting for 
the associated systematic uncertainty, we likewise determine the significance of the evidence 
for the decay B — > x C 2(direct)X to be only 1.4a. 

To calculate the branching fractions B{B — > XcipX), we use the measured signal yields 
N(B — > Xci,2X), the reconstruction efficiencies, the number of produced BB pairs, and the 
daughter branching fractions. The reconstruction efficiencies, determined from simulation, 
are (25.7 ± 0.2)% for Xd and (26.6 ± 0.2)% for x c2 , where the uncertainties are due to the 
size of our B — ► XcipX simulation samples. For the calculation of the rates for the decays 
B — > Xci,2 (direct) X, we make an assumption that the only other source of Xd,2 production 
in B decays is the decay chain B — >■ ip(2S)X with ip{2S) — > Xa,2l- The 95% confidence 



intervals are calculated using the Feldman- Cousins approach ||16|| . The resulting branching 
fractions are listed in Table |. Taking into account correlations between the uncertainties, we 
obtain the branching ratio T[B -► Xc2(direct)X]/r[ J B -> Xci(direct)X] = 0.18 ± 0.13 ± 0.04; 
the 95% CL upper limit on the ratio is 0.44. 



TABLE I. Branching fractions for inclusive B decays to Xcl an d Xc2- 

Branching fraction Measured value 95% CL interval 

(xl0~ 3 ) (xlO" 3 ) 

B(B -> XdX) 4.14 ±0.31 ±0.40 — 

B[B ^> Xd(direct)X] 3.83 ± 0.31 ± 0.40 

B{B->Xc2X) 0.98 ±0.48 ±0.15 [0.2,2.0] 

B[B -* x C 2(direct)X] 0.71 ± 0.48 ± 0.16 [0.0,1.7] 



The systematic uncertainties are listed in Table [D]. The sources of the uncertainty can 
be grouped into three categories: 

Fit procedure. — This category includes the uncertainties due to our choice of the signal 
and background shapes as well as the bin size. To fit the x& an d Xc2 signal, we use the 
templates extracted from simulation. We therefore are sensitive to imperfections in the 
simulation of the photon energy measurement. The systematic uncertainties associated with 
the simulation of the calorimeter response are estimated by comparing the 7r° — > 77 invariant 
mass lineshapes for inclusive 7r° candidates in the data and in Monte Carlo samples. Then 
the Xci and Xc2 templates are modified accordingly in order to determine the resulting 
uncertainty in the signal yields. To estimate the uncertainty associated with the calorimeter 
energy scale, we shift the Xci and Xc2 templates by ±0.6 MeV/c 2 in the fit. The uncertainty 
due to time-dependent variations of the calorimeter energy scale is small compared to the 
overall energy scale uncertainty. To estimate the uncertainty due to the calorimeter energy 
resolution, we change the width of the xa and Xc2 templates by ±4%. The uncertainty in 
the background shape is probed by fitting the background with a template extracted from 
high-statistics samples of simulated T(45) — ► BB and non-55 continuum events; only the 
template normalization, not its shape, is allowed to float in the fit. 

Efficiency calculation. — This category includes the uncertainties in the number of pro- 
duced BB pairs, tracking efficiency, photon detection efficiency, lepton detection efficiency, 
and model-dependence and statistical uncertainty of the B — > XcipX simulation. The Xc\.2 
polarization affects the photon energy spectrum. We define the helicity angle Oh to be the 
angle between the 7 direction in Xc frame and the Xc direction in the B frame. We assume 
a flat cos 9h distribution in our simulation. The systematic uncertainty associated with this 
assumption is estimated by comparing the reconstruction efficiencies in the Monte Carlo 
samples with I (Oh) oc sin 2 9h and I (Oh) oc cos 2 Oh angular distributions. Parity is conserved 
in the decays Xci,2 ~> J l^li so the helicity angle distribution contains only even powers of 
cosOh- Another source of uncertainty is our modeling of the X system in the B — > Xci^X 
simulation. Photon detection efficiency depends on the assumed model through the Xc mo- 
mentum spectrum and the 7r° multiplicity of the final state. In our simulation, we assume 
that X is either a single K or one of the higher K resonances; we also include the decay chain 



B —* ip(2S)X with ip(2S) — ► Xci,2l- To estimate the systematic uncertainty, we compare the 
Xc -^ J/^l detection efficiency extracted using this sample with the efficiency in the sample 
where we assume that X is either a K ± or Kg — > 7r + 7r~. 

Assumed branching fractions. — This category includes the uncertainties on the exter- 
nal branching fractions. We use the following values of the daughter branching fractions: 
B(J/i/j -> £ + £~) = (5.894 ± 0.086)% 0, B( X d -> J/^l) = (27.3 ± 1.6)% 0, and 



B{Xc2 -> J/Vt) = (13.5 ± 1.1)% 0. In the calculation of B[S -> x c i )2 (direct)X], we 
also assume the following values: B(B -► i/j(2S)X) = (3.5 ± 0.5) x 10" 3 0, B(ip(2S) 



Xdl) = (8.7 ± 0.8)% 0, and B(if>(2S) -> x c27 ) = (7-8 ± 0.8)% g 



TABLE II. Systematic uncertainties on B(B — > Xci,2X). 



Source of relative uncertainty in % 

systematic uncertainty B(B — > Xcl-^O £>(!? — ► Xc2-^) 

Fit procedure 

7 energy scale 0.4 5.6 

7 energy resolution 2.8 6.9 

Background shape 1.8 6.8 

Bin size 0.0 1.9 

Efficiency calculation 

N{BB) 2.0 2.0 

Tracking efficiency 2.0 2.0 

Lepton identification 4.2 4.2 

Photon finding 2.5 2.5 

Monte Carlo statistics 0.7 0.7 

Model for linB^ Xci,2* 3.3 3.3 

Polarization of Xd,2 1-0 1.0 

Assumed branching fractions 

B(Xd,2 -» J I in) 5.9 8.1 

B(J/ip -► r^-j 1.5 1.5 

B(S -» V(2^)^) a 1.1 5.5 

g(^(25) -> Xd,2lT 07 4.0 

a Contributes only to uncertainty on B[B — > Xcr,2(dfrect)X]. 

In the second part of this work, called the 5-reconstruction analysis, we employ the 5- 
reconstruction technique similar to the one developed for the b — > 57 rate measurement [[17] . 
We still extract Xci and Xc2 signal yields from a fit to M( J/^7) — M(J/ip) distribution, 
but we select only those J/ip'y combinations that reconstruct to a B — > J/ip^fX a decay. 
This ^-reconstruction technique is used to suppress backgrounds and allows us to probe 
the composition of the X s system accompanying Xc\,2 mesons. We extract the branching 
ratio Ti{Xc2/Xci) = F[B — > x c2 (direct)X s ]/r[i? — > x cl (direct)X s ] for the following three X s 
configurations: 

1. Sample A. — X s is reconstructed as a kaon (K + or K s — > 7r + 7r^) with to 4 pions, 
one of which can be a tt°. We consider 21 possible X s modes as well as the charge 
conjugates of these modes. 

7 



2. Sample B. — X s is reconstructed as a single kaon or if* (892). A Kit combination is 
a K* candidate if \M(Kn) — Mk*\ < 75 MeV/c 2 , where Mk* is the world average 
K*(892) mass fTJ. 

3. Sample C. — X s is reconstructed as a kaon with 1 to 4 pions, but not as a i^*(892) 
candidate (\M(Kit) - M K *\ > 200 MeV/c 2 ). 

Thus samples B and C are subsets of A. To an excellent approximation, sample A is a 
sum of B and C. With sample A, we try to reconstruct as many B — > J/ip^Xg decays as 
possible. Dividing sample A into subsamples -B and C, we also probe the dynamics of the 
B — > Xci,2^s decays. If the dominant production mechanisms for Xci and Xc2 are different, 
color-singlet mechanism for Xci and color-octet for x c2 , then it is natural to expect that Xc2, 
in comparison with Xci, is more often accompanied by multi-body X s states rather than a 
single K or K*. Thus the measured Xc2~to-x c i production ratio might be quite different for 
samples B and C. 

We require that the charged kaon and pion candidates have, if available, dE/dx and time- 
of-fiight measurements that lie within 3a of the expected values. The dE/dx measurement 
is required for kaons, but used only if available for pions. The time-of-flight measurement 
is used only if available. The K° s — > 7r + 7r~ candidates are selected from pairs of tracks 
forming displaced vertices. We require the absolute value of the normalized Kg — > tt + 7t~ 
mass to be less than 4 and perform a fit constraining the mass of each K s candidate to 
the world average value |14]]. Photon candidates for ir° — > 77 decays are required to have 



an energy of at least 30 MeV in the central region and at least 50 MeV in the endcap 
region (0.71 < |cos6> 7 | < 0.95) of the calorimeter. We require the absolute value of the 
normalized 7r° — > 77 mass to be less than 3 and perform a fit constraining the mass of each 
7T° candidate to the world average value |14]. The J/i/j four- momentum used in B — > J/tp^yXs 



reconstruction is obtained by performing a fit constraining the J/ip candidate mass to the 
world average value ]rj[]. 

The B candidates are selected by means of two observables. The first observable is the 
difference between the energy of the B candidate and the beam energy, AE = E(B) — -Bbeam- 
The average AE resolution varies from 12 to 17 MeV depending on the ^-reconstruction 



mode. The second observable is the beam-constrained B mass, M(B) = \JE^ eam —p 2 (B), 
where p(B) is the B candidate momentum. The average M(B) resolution is 2.7 MeV/c 2 and 
is dominated by the beam energy spread. We use the normalized M(B) and AE variables 
and require \AE\/a(AE) < 3 and \M(B) - M B \/a(M) < 3, where M B is the nominal B 
meson mass. The fit to M(J/tp / -f) — M{J/ip) distribution is then performed in the same 
manner as in the inclusive analysis. We still use a 5th order Chebyshev polynomial to fit 
the background for samples A and C, but we reduce the order of the polynomial to 3 for the 
low-statistics sample B. The fits are shown in Fig. [J] and the Xci and Xc2 signal yields are 
listed in Table |T|. The S-reconstruction technique renders negligible the contribution from 
n.on.-BB continuum events. We finally subtract the ip{2S) — > Xci,27 feeddown to obtain the 
rates for direct Xa,2 production in B decays. For all three X s configurations, we observe a 
strong Xci signal yet no statistically significant signal for direct Xc2 production (Table fl|). 
To calculate the branching ratio 7l(Xc2/Xci), we multiply the ratio of the feeddown-corrected 
Xci, 2 yields by the reconstruction efficiency ratio S(x c i)/^(Xcn) and by the branching ratio 



r(Xci ~~ * J/' l Pl)/T(Xc2 - ► J/ifrl)' The efficiency of the B- reconstruction depends on the 
composition of the X s system. We assume that the X s system composition is the same for 
Xd and Xc2 production. From our simulation we determine £{Xci)/£{Xca) — 0-93 for all 



three X s configurations. The resulting Xc2-to-x c i production ratios are listed in Table III. 



TABLE III. Results for each of the three X s configurations used in B — ► J/ip^yXg reconstruc- 
tion. The Xcl an d Xc2 event yields with associated statistical uncertainties are listed in lines 1 and 
2. Line 3 contains the significance of the B — > Xc2(direct)X s signal with statistical and systematic 
uncertainties taken into account. Lines 4 and 5 contain the measured value and 95% confidence 
interval for the branching ratio 1Z(Xc2/Xd) = T[B — > Xc2(diiect)X s ]/T[B — > Xci(dfrect)X s ], deter- 
mined with an assumption that the X s system composition is the same for x& and Xc2 production. 





Sample A 


Sample B 


Sample C 


N(B -» XciX s ) 


279 ± 25 


96 ±12 


183 ± 22 


N(B -» Xc2X s ) 


oi+rs 


13-9±£S 


18 ±16 


Significance of B — > Xc2(direct)X s 


1.2a 


2.0a 


0.6a 


KiXa/Xci) 


0.18 ±0.12 ±0.09 


0.27±{j:i3 ± 0.05 


0.14 ±0.18 ±0.14 


95% CL interval for Il{xc2/Xci) 


[0.00,0.48] 


[0.04,0.58] 


[0.00,0.59] 



The systematic uncertainties for the ^-reconstruction analysis are listed in Table 0. 
The sources of uncertainty can be grouped into the following four categories: 

Fit procedure. — As in the inclusive analysis, we estimate the uncertainties in the signal 
and background shapes. We shift the Xa,2 templates by ±0.6 MeV/c 2 and vary their widths 
by ±4%. The requirement on AE in B — > J/ifj^Xg reconstruction truncates the low-side 
tail of the Xci,2 shapes. We estimate the uncertainty due to this effect by using the Xc\.2 
templates obtained from the simulation with a requirement that the measured Xc energy is 
within 3<T of the generated value. The uncertainty in the background shape dominates the 
fit procedure uncertainty. To probe this uncertainty, we fit the background with different 
templates, allowing only the template normalization, not its shape, to float in the fit. One 
template is extracted from simulation separately for each of the samples A, B, and C. 
Another template, the same for all three X s configurations, is the background shape from 
the inclusive analysis (Fig. |l|a). 

'0(25') subtraction. — - The sources of the systematic uncertainty associated with the 
■0(2S')-feeddown subtraction include B{B — >■ if;(2S)X), B{jp{2S) — > Xci,2j), the size of our 
B — > ip(2S)X simulation sample, and the composition of X in B — > ip(2S)X decays. To 
estimate the uncertainty due to our model of the X system composition in the B — > ip(2S)X 
simulation, we check whether the data and the simulation agree on the ratio of ip{2S) — * £ + £~ 
event yields obtained in the inclusive reconstruction and after the B — ► tp(2S)X s reconstruc- 
tion. This category also includes the uncertainties that would have canceled for the ratio 
TZ(xc2/Xci) were it not for the V ; (25')-feeddown subtraction. These sources of uncertainty 
are B{J/ip — ► £ + £~), N(BB), tracking, photon finding, and lepton identification. 

£(Xc2)/£(xa)- — We assume that the X s system in B — > Xa,2X s is the same for xa 
and Xc2- We do not assign any uncertainty for this assumption. The remaining sources 
of uncertainty are the Xci,2 polarization and the statistics of the B -^ Xa,2X s simulation 
samples. 
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B(Xd,2 —* J/i>l)- — Our measurement depends on the ratio r(% cl —* J/ipj)/T(x C 2 —* 
J/if>j) and its uncertainty. 

TABLE IV. The absolute systematic uncertainties on the branching ratio T^iXcu/Xcl) f° r each 
of the three X s configurations used in B — > J/ip^fX s reconstruction. 

uncertainty on Tl(xc2/Xci) 

Source of uncertainty Sample A Sample B Sample C 

Fit procedure 0.084 0.039 0.142 

ip(2S) subtraction 0.007 0.001 0.006 

£{Xci)/£{Xc2) 0.003 0.006 0.003 

g(Xci,2 -> J/ih) 0.022 0.026 0.019 

Added in quadrature 0.09 0.05 0.14 



In conclusion, we have measured the branching fractions for inclusive B decays to the 
Xci and Xc2 charmonia states. Our measurements are consistent with and supersede the 
previous CLEO results ||. We have also studied B — > XcipXs decays, reconstructing 
X s as a kaon and up to four pions. In this way, we have measured the branching ratio 
T[B — > Xc2(diiect) X S ]/T[B — > Xci(diiect)X s ] for three X s configurations. In all the cases, 
we observe strong xa signal yet no statistically significant signal for Xc2 production. Our 
measurement of the Xc2-to-x c i production ratio in B decays is consistent with the predic- 
tion of the color-singlet model [0] and disagrees with the color-evaporation model [EJ. In 
the NRQCD framework, our measurement suggests that the color-octet mechanism does not 
dominate in B — > XcX decays. 
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FIG. 1. The M^J/ipj) — M(J/ip) distribution in the Y(4S) data (points with error bars). 

Plot (a) is for inclusive J/ip^f combinations, whereas plots (b), (c), and (d) are for those J/tp'f 

combinations that reconstruct to a B — > J /ip^X s decay with the X s composition corresponding 

to samples A, B, and C described in the text. The fit function is shown by a solid line with the 

background component represented by a dashed line. The insets show the background-subtracted 

distributions with the Xd and Xc2 fit components represented by a solid line. 
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